Crossover from Goldstone to critical fluctuations: Casimir forces in confined O(n) 

symmetric systems 
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We study the crossover between thermodynamic Casimir forces arising from long-range fluctua- 
tions due to Goldstone modes and those arising from critical fluctuations. Both types of forces exist 
in the low-temperature phase of 0(n) symmetric systems for n > 1 in a d-dimensional L^' 1 x L 
slab geometry with a finite aspect ratio p = L/Lu. Our finite-size renormalization-group treatment 
for periodic boundary conditions describes the entire crossover from the Goldstone regime with a 
nonvanishing constant tail of the finite-size scaling function far below T c up to the region far above 
T c including the critical regime with a minimum of the scaling function slightly below T c . Our 
analytic result for p <C 1 agrees well with Monte Carlo data for the three-dimensional XY model. 
A quantitative prediction is given for the crossover of systems in the Heisenberg universality class. 

PACS numbers: 05.70.Jk, 64.60.-i, 75.40.-s 



In the past two decades, substantial effort has been de- 
voted to the study of thermodynamic Casimir forces [l| 
that result from two fundamentally different sources in 
confined condensed matter systems: (i) from classical 
fluctuations with long-ranged correlations due to mass- 
less Goldstone modes [2|, |3J, and (ii) from long-ranged 
critical fluctuations Both types of fluctuations exist 
in the low-temperature phase of O(n) symmetric systems 
such as superfluids @, superconductors [f|, XY mag- 
nets (n = 2) @, 0, @|, and isotropic Heisenberg magnets 
(n = 3) Q. While successful analytic theories have been 
developed that separately describe such Casimir forces ei- 
ther (i) in the Goldstone-dominated regime deeply in the 
low-temperature phase (2, 3] or (ii) above bulk criticality 
[9l4T^| , there is a serious lack of knowledge concerning an 
analytic theory of the crossover between these two types 
of Casimir forces in the low-temperature phase. The goal 
of this Letter is to provide such a crossover theory for the 
case where the confining geometry has periodic boundary 
conditions. Such systems are well accessible to numeri- 
cal studies 0, 0, We shall present analytic results 
for general n that are in very good agreement with the 
existing Monte Carlo (MC) data for n = 2 both in the 
Goldstone regime as well as in the critical region includ- 
ing the crossover between these regions. We also present 
a quantitative prediction for this crossover of the Casimir 
force in the (n — 3) Heisenberg universality class. The 
concept of our theory should be applicable also to the 
case of Dirichlct boundary conditions which are relevant 
to the crossover of Casimir forces in superfluids @ and 
superconductors 0]. 

A crucial ingredient of our approach for general n is an 
appropriate choice of the geometry. We consider a finite 
d dimensional x L slab geometry with a finite aspect 
ratio p = L/L\\. This is well justified by the fact that all 
experiments and computer simulations were performed 
in slab geometries with a small but finite p rather than 
p = 0. Recent Monte Carlo data for p — 1/6 [3j and 
p = 0.01 @ show that the p - dependence of the Casimir 
force is quite weak for p -c 1. The basic advantage of a 
finite-slab geometry (rather than the idealized oo^ 1 x L 



film geometry studied earlier H, H, ) 1S * ne ab- 

sence of singularities of the free energy density at finite 
temperatures and the existence of a discrete mode spec- 
trum with a dominant lowest mode that is amenable to a 
simultaneous analytic treatment of the low-temperature 
and the critical regions. 

Our analytic treatment is based on the O(n) symmetric 
isotropic (p Hamiltonian 



(1) 



where </?(x) = V~ 1 £kV3ke lkx is an n-component field in 
a finite volume V = L'^~ 1 L with periodic boundary con- 
ditions. The summation Ek runs over discrete k vectors 
including k = up to some cut off A. The fundamen- 
tal quantity from which the Casimir force per unit area 
Fcas = — d[Lf cx ]/dL can be derived is the excess free 
energy density (divided by UbT) f cx — f — fb where 



f(T,L,L« 



V -1 ln / Vipexp(~H) 



(2) 



and fb = limy^oo / are the free energy densities of the 
finite system and the bulk system, respectively. 

It is expected that, for isotropic systems near criticality 
and for large L and L», Fq (is can be written in a finite- 
size scaling form [l5| 



F Cas (i,L,L||) = L- d X(i,p) 



(3) 



with the scaling variable x = t{L/^Q + ) 1 / v , t = (T — 
T c )/T c where £o+ is the amplitude of the bulk correla- 
tion length above T c . As noted for the case of film ge- 
ometry (p — > 0) [3], the scaling form ^ applies to the 
low-temperature region where the scaling function satu- 
rates at a nonzero negative value X(— oo, 0) < for n > 1 
[US EH- S° far, however, no analytic calculation of the 
function X(x,0) for systems with n > 1 and periodic 
boundary conditions has been performed that describes 
X in the whole low-temperature region — oo < x < [l6| . 
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It is the goal of this paper to describe the full crossover 
in terms of a single scaling function X(x,p), for general 
n > 1 and small p, from the Goldstone-dominated be- 
havior for x — > — oo up to the high-temperature behavior 
for J ^ 1 including the critical region < |x| ~ 0(1) 
above and below T c . 

An important conceptual difference between the previous 
perturbation approach within the ip theory for {n = 1) 
systems with a discrete (Ising-like) symmetry [ll[ and for 
systems with a continuous symmetry (n > 1) presented 
in this paper is the following. For n = 1, two sepa- 
rate perturbation approaches were necessary for describ- 
ing finite-size effects both in the central finite-size regime 
near criticality and those in the low-temperature phase 
far below T c in order to capture the two-fold (spin up and 
spin down) degeneracy characteristic of the ground state 
of Ising-like systems. For n > 1, there exist both longitu- 
dinal and transverse fluctuations. The latter correspond 
to orientational changes that permit the order parameter 
fluctuations to exhaust the full phase space since there 
exists no large barrier like that between the purely lon- 
gitudinal spin-up and spin down configurations of large 
Ising-like systems far below T c . As a consequence, a sin- 
gle perturbation ansatz suffices for n > 1 to capture both 
the critical fluctuations and the fluctuations due to the 
Goldstone modes. This results in a smooth description 
of the crossover between the two different regimes with- 
out the necessity of matching two separate pieces of the 
theory. 

We decompose tp = $ + a into a homogeneous lowest- 
mode $ = V -1 J d xcp(x) and higher-mode fluctua- 
tions er(x) = V r ~ 1 Ek#o<^ke lkx . We further decompose 
cr(x) into "longitudinal" and "transverse" parts cr(x) = 
<tl (x) + ot (x) which are parallel and perpendicular with 
respect to $. Correspondingly, the Hamiltonian H is 
decomposed as H = H (&) + H(&, a) where 



H m = v 



i 



-r $ 2 + u ($ 2 )' 



The second approximation is to neglect the effect of the 
higher-mode parts oc Qa^a 2 and cx (c 2 ) 2 on the finite- 
size properties (but not on the bulk critical exponents 
which will be incorporated via Borel-resummed field- 
theoretic functions). No further approximation (such 
as an e = 4 — d expansion) will be made in the sub- 
sequent renormalization-group treatment at fixed d. Our 
approach goes beyond a naive Gaussian approximation 
not only because of the fourth-order term ~ uq{<& 2 ) 2 
in Hq but also because of the size- dependent couplings 
~ uoMqCTl 2 and ~ uqM 2 ut 2 between the lowest-mode 
average and the higher modes which arise from the terms 
^olOl 2 and r~oTOT 2 , respectively. After integration over 
<7, we obtain the unrenormalized free energy density 



/ = / - i In { j d n $ exp [-#„(*)] } 



(6) 



with So(r) = V 1 Sk^o m ( r + k 2 ) where /q is indepen- 
dent of tq and uq. 

For finite V, the transverse parameter foT remains posi- 
tive in the whole range — oo < ro < oo and interpolates 
smoothly between the large- volume limit above and be- 
low T r , 



(4) 



lim roT = 



ro 




for 
for 



r > 
r < 



(7) 



H(*,*) = J d d x{ i [Volol 2 + roT^T 2 + (Vol) 2 
v 

+ (Vcr T ) 2 ] + 4?i $criV + u Q (a 2 ) 2 } (5) 

with the longitudinal and transverse parameters 
r L($ 2 ) = r + 12u $ 2 , r 0T ($ 2 ) = r + 4u $ 2 . The free 
energy density / is then calculated by first integrating 
over a and subsequently over $. 

An important reference quantity of our theory 
is the n-dependent lowest-mode average M 2 — 
Jd n $$ 2 exp[- J ff ($)] / Jd n <f>exp[-H ($)]. The phys- 
ical importance of the quantity M5 2 is related to the 
fact that the main contribution of the integration over 
$ comes from the region around <£ 2 « Mq and that this 
quantity is relevant in the whole range — 00 < ro < 00 
above and below T c . This provides the justification 
for replacing the $ 2 -dependent parameters roL( ( & 2 ) and 
r 0T ($ 2 ) in H by f 0L = r 0h (M$) and f 0T = r 0T (Af 2 ). 



The vanishing of the transverse parameter tqt is the char- 
acteristics of the massless Goldstone modes which are the 
origin of long-range correlations and the Casimir force 
well below T c . In higher-order perturbation theory for 
the bulk system, spurious (infrared) singularities arise 
due to the vanishing of tot [171 ]. Within our approxima- 
tion, such spurious singularities do not yet appear since 
S'o(roT) has a finite large- volume limit below T c . 

The bare expression ([6]) does, of course, not yet cor- 
rectly describe the crossover from the Goldstone to the 
critical regime. Both additive and multiplicative renor- 
malizations are necessary, after subtracting a nonsingular 
bulk part f ns . Integration of the renormalization-group 
equation then leads to the scaling form f s (t,L,Lu) = 
L~ d F(x,p) for the singular part f s = f — fns- We have 
performed these steps within the minimal renormaliza- 
tion scheme at fixed dimensions 2 < d < 4 flal. The 
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result reads 



F(x,p) = 


= - An 


4d " 


(n-1) 


A Te 


]d/2~ 

~d~ 


+ /"'{ 


-5- 






-m(2 J 




1 exp 










2o 



B(u*) 



oo ^ 

J (z,p) = p 1 - d J dyy-H exp [-z2,/(4tt 2 )] 



(8) 



(ir/yf/ 2 - pK{p 2 y) K(y) + p 



-p d 1 exp(-y) 



„d-l 



(9) 



with if (y) = X)m=-oo exp(-ym 2 ) and with 

h(2,p) = P 8i^u^ 2 p( d -^A- d 1/2 My), (io) 



(11) 



^=jds^^,jds^*-°- 

o o 
where l(x,p) and y(x,p) are determined implicitly by 

y + 120 2 (y) = ^"^/^V V" 1 / 2 , (12a) 
y = 5 Q* |- a /( 2 ^( 1 - (i )/ 2 4y 2 U , - 1/2 , (12b) 

The fixed point value of the renormalized four-point 
coupling u is denoted by u* . The quantities Q* = 
Q(l,u*,d) and B(u*) are the fixed point values of the 
n dependent amplitude function Q(l, u, d) of the second- 
moment bulk correlation length above T c and of the field- 
theoretic function B{u) related to the additive renormal- 
ization of bulk theory [3], respectively. Ad = r(3 — 
d/2)[2 d ~ 2 7r d / 2 (d — 2)] _1 is an appropriate geometric fac- 
tor. 

Eq. ijHJ) is valid for general n and 2 < d < 4 above, at, 
and below T c including the Goldstone regime for n > 1. 
It incorporates the correct bulk critical exponents a and 
f and the complete bulk function B(u*) (not only in one- 
loop order). F(x, p) is an analytic function of x at finite p, 
in agreement with general analyticity requirements. For 
n = 1, the function © is not applicable to the region well 
below T c as it does not capture the two-fold degeneracy 
of the ground state of Ising-like systems, as discussed in 

n. 
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FIG. 1: (Color online) Scaling function X(x, p) of the Casimir 
force as a function of x — t(Z//£o+) 1 ^ 1 ' (a) for n = 2 and (b) 
for n = 3 in three dimensions. Solid lines: d — 3 theory for 
p — 1/6 as calculated from ((8)). (|13|) . Dot-dashed and dotted 
lines: e expansion results for p = 0, Ej. Horizontal dashed 
lines: Casimir amplitude X( — oo,0) = — (n — 1)((3)/tt. MC 
data (a) from [3j for the XY model and (b) from [7J for the 
Heisenberg model. 



To obtain the scaling function X of the Casimir force we 
consider the singular part of the bulk free energy density 
ff b (f) = A ± \t\ dl '. It can be written as f± b = L' d F b ± {x) 

where F b (x) is the bulk part of F(x,p). It is derived 
from ||5J| in the limit of large \x\. The scaling function of 
f cx is then given by F cx (x,p) — F(x,p) — F, ,{x) from 
which we calculate the desired scaling function 

x(x, P) = (d-i )F -(x,p) - ^j^ia -ftqmu) 

v ox op 

In the limit p —¥ 0, x — > — oo, our function X(x, p) yields 
the finite low-temperature amplitude 

X{-oo, 0) = -(n - l)(d - l)ir- d/2 r{d/2)((d) (14) 

of the Goldstone regime in film geometry. For n = 2,d = 
3 this agrees with Eq. (25) of Q. 

MC data are available for both the three-dimensional 
XY [3j and Heisenberg [7] models with the aspect ratio 
p = 1/6. The comparison of our result for p = 1/6 with 
the MC data is shown in Figs. 1 (a) and (b) for n = 2 
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and n = 3, respectively. Here we have employed the fol- 
lowing numerical values [nj [2(j: v = 0.671,0.705, a — 
-0.013,-0.115, u* = 0.0362,0.0327, Q* = 0.939,0.937, 
B(u*) = 1.005,1.508, for n = 2,3, respectively. Also 
shown are the e expansion results for p = from @, QjJ 
for T >T C (dot-dashed and dotted lines). The horizon- 
tal dashed lines represent the p = Casimir amplitude 
(|14p due to the Goldstone modes. For n = 2 there is 
an excellent overall agreement of our result (solid line) 
with the MC data Q including the crossover from the 
Goldstonc-dominated region to the critical region. 

For n — 3 [Fig. 1 (b)] there are systematic deviations 
between all theoretical curves and the MC data of [7j . A 
similar discrepancy with the MC data of Q exists also 
for n = 2 [not shown in our Fig. 1 (a)l as shown in Fig. 6 
of [H for the XY model. As noted in Q , these discrepan- 
cies might be due to the uncertainty in the normalization 
factor used for the MC data in [7J. More accurate MC 

fo: 
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FIG. 2: (Color online) Critical amplitude F(Q,p), ©, for 
7i = 2 and d — 3 at T c (solid line). The maximum —0.3275 is 
at pmax = 0.196. The dashed line is the extrapolation from 
p = Pmax to p — (film geometry). The dotted line represents 
© at T c for p < pmax ■ MC estimate (triangles) for the d — 3 
XY model for p = 0.01, 1/4, 1/2 from ^ . 

Finally we discuss the question as to the expected range 



of applicability of our theory at finite p. For this purpose 
we consider the finite-size amplitude F(0, p) of Eq. © at 
T = T c . We expect a monotonic p - dependence of F(0, p) 
on the basis of the monotonicity hypothesis inferred from 
previous results for n = 1 [HI, HH and n = oo [ll| . A 
comparison with MC data for the three-dimensional XY 
model [22] is shown in Fig. 2 in the range < p < 1/2. 
The MC data show a very weak p - dependence for p <C 1 . 
Our result is in reasonable agreement with the MC data, 
except that the monotonicity of the p - dependence is not 
reproduced by our theory for small p < p max = 0.196 
(dotted portion of the curve). The deterioration of the 
quality of our theory for p — > is to be expected since 
the separation between the lowest mode and the higher 
modes goes to zero in this limit. For this reason it is 
expected that Eqs. ([5]) - (fT2"]) are not quantitatively re- 
liable for p <C 0.2 in some region \x\ -c 1 close to T c . 
In particular, the scaling function of the film free energy 
density Ffu m (x) — limp^o F(x, p) obtained from ([5]) is 
not expected to be reliable for \x\ -C 1. This function 
has, in fact, an artificial cusp-like singularity at x — 
similar to that of previous approximate p — theories 
@, [t| [13, fl3 - IT3 |. Our monotonicity criterion provides 
a theory-internal argument against the reliability of ap- 
proximate p — results in the region \x\ <C 1. Neverthe- 
less, assuming a negligible p - dependence for p < p max 
we obtain from the extrapolation of the maximum in Fig. 
2 to p = (dashed line in Fig. 2) our prediction 



F iUm (0) « F{0 

i Pmax ) 



-0.3275 



(15) 



for film geometry at bulk T c . This is in acceptable agree- 
ment with the MC estimate -0.292 at p = 0.01 [22] 
shown in Fig. 2. 

It would be interesting and important to extend our ap- 
proach to the most relevant case of Dirichlet boundary 
conditions in order to explain the experimental results 
for the Casimir force in superfluid 4 He Q for which there 
exists no satisfactory analytic theory so far. 

The author is grateful to D. Dantchev, M. Hasenbusch, 
and O. Vasilyev for providing the MC data of % 0, HI 
in numerical form. 
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